GROUND STATE PROPERTIES OF GRAPHENE IN 
HARTREE-FOCK THEORY 

CHRISTIAN HAINZL, MATHIEU LEWIN, AND CHRISTOF SPARSER 

Abstract. We study the Hartree-Fock approximation of graphene in infi- 
nite volume, with instantaneous Coulomb interactions. First we construct its 
translation-invariant ground state and we recover the well-known fact that, due 
^~v ' to the exchange term, the effective Fermi velocity is logarithmically divergent 

»vj , at zero momentum. In a second step we prove the existence of a ground state 

in the presence of local defects and we discuss some properties of the linear 
response to an external electric field. All our results are non perturbative. 
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Graphene is a mono-crystalline graphitic film in which electrons behave like 2- 
j^ , dimensional Dirac fermions without mass. It has attracted a huge interest in the 

last decades [H] . The Fermi surface of metals in dimension d is usually (the union 
of) {d— l)-dimensional manifold(s), whereas for gapless semi-conductors it is often 
composed of parabolic-like points. Graphene, on the other hand, has a peculiar 
quasi-particle dispersion relation which is conical at the degeneracy points, leading 
to the effective massless Dirac equation. This is not so exceptional, however. It was 
recently shown that conical singularities are generic in (non-relativistic) quantum 
crystals having the honeycomb lattice symmetry [13]. The conical dispersion rela- 
tion, combined with the fact that the Fermi velocity v-p is 300 times smaller than 
the speed of light, makes graphene an ideal condensed matter system for testing in 
the lab our understanding of 2D massless relativistic particles [42l |48] . 

Quantum Electrodynamics (QED) is a very powerful theory which, however, 
has only been rigorously formulated in a perturbative fashion. The small value 
of the Fermi velocity in graphene restricts the validity of perturbation theory and 

1 



2 C. HAINZL, M. LEWIN, AND C. SPARBER 

it is important to be able to resort to non-perturbative methods. Short range 
interactions in graphene are well understood [3T1 [TCI [H] . They do not modify the 
general properties of the ground state as compared to the free case. The effect 
of long range Coulomb interactions, however, is much less clear. It was argued 
in [301110] that instantaneous Coulomb interactions lead to an additional logarithmic 
divergence in the infrared regime for the interacting Fermi velocity, whereas the 
authors of [H] claimed that retarded effects make it tend to the speed of light. 
In [TS] a careful renormalization group analysis for retarded interactions indicated 
that the Fermi velocity actually has a limit, but that it is in general smaller than 
the speed of light. 

In this work we shall not consider retardation effects, but instead confirm the 
logarithmic divergence in the case of instantaneous Coulomb interactions, in a fully 
non-perturbative setting, based on a Hartree-Fock model. In [251 1261 129 ) . two of 
the authors of the present paper have, with E. Sere and J. P. Solovej, studied the 
Hartree-Fock approximation of 3-dimensional Quantum Electrodynamics. These 
works were inspired by an important paper of Chaix and Iracane ^6 and they have 
been summarized in |28[ 112) . They are purely non-perturbative and hold for all 
values of the coupling constant < a < i/w ~ 1.27. For a > 4/7r, the system 
is known to become unstable [7]. In the present work we extend these results to 
(massless) electrons in two dimensions. We identify the exact Hartree-Fock ground 
state of the system at half filling and zero temperature under the sole assumption 
that 

< a < 0.48637 or, equivalently, wp > 2.0560, 

which covers the case of graphene (see Section [2] below). 

Our methodology is as follows. We consider a Hartree-Fock type model in which 
particles interact through the instantaneous Coulomb potential and with a kinetic 
energy given by the massless Dirac operator. Since we do not use normal-ordering, 
the Hamiltonian is unbounded from below. However, we shall, as a first step, 
construct the free Dirac sea (i.e., the absolute minimizer of the energy in the absence 
of external fields) by means of a thermodynamic limit. This state corresponds to 
filling the negative energies of an effective mean-field translation-invariant operator 
of the form Vesip)^'' • P- Here, Vesip) denotes an effective Fermi velocity which we 
shall compute exactly and show that, for small momentum, it diverges like 

vcs{p) -^ tIoSm' 
P-i-o 4 \p\ 

where A is a fixed ultraviolet cut-off. In a second step, we shall introduce an external 
electrostatic field and obtain a bounded-below energy by subtracting the (infinite) 
energy of the free Dirac sea. This enables us to prove the existence of a ground 
state in the presence of the external field, in infinite volume. In other words, we 
use the translation-invariant free state as a reference and we describe variations 
compared to it. 

In addition, we shall calculate the linear response function B{p) of the density 
of charge of graphene to an external electrostatic field, and we shall show that, in 
the limit of small momenta, 

B{p) 



'^° 4 log (^ 

It is interesting to note that the dielectric behavior of Hartree-Fock graphene is 
universal at low momentum, that is, independent of the bare Fermi velocity vp. As 
we shall see, this follows from the fact that the effective Fermi velocity Ves(,p) has 
the same property. 
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As compared to the previously quoted results in 3 dimensions, we deal here with 
massless Dirac fermions. Mathematically, this creates a certain lack of control of 
the behavior of minimizing sequences in the infrared domain, which complicates the 
study of the exchange term. Interestingly, though, our existence proof for ground 
states heavily relies on the effective Dirac operator including the logarithmically 
effective Fermi velocity, and which can be used to get a better control. On the 
other hand, the ultraviolet behavior which was problematic in the three dimensional 
case is not an issue here. For graphene it would even have no meaning since the 
Dirac dispersion relation is only valid for low momentum anyhow. In this paper we 
consequently impose a sharp ultraviolet cut-off A which is kept fixed all along our 
study and which mimics the presence of the underlying lattice. 

The paper is organized as follows. In the next section we properly introduce the 
Hartree-Fock approximation of massless 2D QED. Then, in Section |3] we construct 
the free Fermi sea of graphene and we discuss the Fermi velocity divergence at 
low momentum due to the Coulomb exchange term. In Section S] we prove the 
existence of a ground state in the presence of external fields, like those induced by 
defects in graphene. This will allow us to compute the linear response of graphene 
in Section [5] Section [5] contains the proof of our main theorem, whereas in the 
Appendix we prove a useful localization formula for the massless pseudo-relativistic 
kinetic operator \/—A. 

Acknowledgment. M.L. acknowledges financial support from the French Ministry 
of Research (ANR-10-BLAN-OIOI) and from the European Research Council under 
the European Community's Seventh Framework Programme (FP7/2007-2013 Grant 
Agreement MNIQS 258023). 

2. HARTREE-FOCK theory of GRAPHENE 

In this section we explain how to derive our model from 2D Quantum Electrody- 
namics, following |29]. We consider a system of 2D massless electrons interacting 
through the 3D instantaneous Coulomb potential. The latter is justified by the 
fact that while the electrons in graphene are essentially confined in 2D, the electric 
field clearly still acts in all three spatial dimensions. The corresponding formal 2D 
QED-type Hamiltonian, written in Coulomb gauge, then reads [501 IHl HI HI] 

(1) H^ = VF / ^*{x) (T ■ (-iV)*(.T) dx + / V{x)p{x) dx 

1 ff p^^)p^y)dxdy 



2JVxr2 \x-y\ 

with <T = ((T^,CT^), the first two Pauli-matrices. Here, and in the following, we shall 
use the notation 

L»" = -ia- ■ V = -ia^d^^ - ia^d^^ 

for the massless 2D Dirac operator. 

We are working in atomic units with the mass and charge of the electrons nor- 
malized to 1. The (bare) Fermi velocity vf is about vp ~ 2.20 in graphene. Thus 
the electrostatic interaction, i.e., the last term of ([T]), and the kinetic part are of the 
same order. The internal variable of the electrons is a pseudo-spin which describes 
to which of the two interlaced triangular lattices of carbon which make up the 
hexagonal lattice the electron belongs. Since we will not consider magnetic fields in 
this paper, we have neglected the (regular) spin of the electron for simplicity. The 
previous Hamiltonian describes well the electrons with energies about the Fermi 
energy of graphene, provided that they live on a much larger scale than the lattice 



4 C. HAINZL, M. LEWIN, AND C. SPARBER 

size. In this regime graphene can be seen as a continuous medium, leading to the 
2D QED-type Hamiltonian ((T|). 

Let us now briefly explain the main objects which enter in the definition of 
the Hamiltonian ([T]). In there ^(x) is the second quantized field operator which 
annihilates an electron at x and satisfies the anti-commutation relation 

(2) **(a;),*(y), + *(2/),**(x), = 2S^^Jix - y), 

with a,i/ G {±1/2} the pseudo-spin variables. The operator p(x) is the density 
operator defined by 

(3) p{x)^2^ ^ , 

cr=l 

where [a, b] = ab — ba. Finally, the function y is a local external electrostatic 
potential which is applied to the system. It is for instance induced by a set of 
defects in the system. The Hamiltonian H^ now formally acts on the fermionic 
Fock space ^ for the electrons. The commutator in the definition ([3]) of p{x) is 
important as it ensures charge conjugation invariance |44) . Precisely, we have 

where "^ is the charge conjugation operator acting on the Fock space. The Hamil- 
tonian MX is unbounded from below on ^ and it is not even a well-defined self- 
adjoint operator. However, it is still possible to define it in a box with suitable 
boundary conditions and with an ultraviolet cut-off, as was done in |29j (see Sec- 
tion [3] below for more details). 

We have already neglected photons in our model. We shall now make another 
approximation, by restricting our attention to Hartree-Fock states. Let us recall 
that the electronic one-body density matrix (two point function) of any electronic 
state fl is defined as 

-f{x,y)^,,, = {^*ix)^'<5j{y),,)n. 
It is an operator on the one-body space such that < 7 < 1, due to the anti- 
commutation relations. It is an orthogonal projection for (and only for) pure 
Hartree-Fock states. In view of ([3]), it is natural to introduce a renormalized one- 
body density matrix 

[*(^);,*(2/).' 



7ron(a;,j/)CT,cr' = ^ri 
By ([2]), we obtain the simple relation 



^ 



I 

7ron = 7-2' 

where / is the identity operator on the one-body space. 

Electronic Hartree-Fock states form a subset of states which are completely de- 
termined by their density matrix 7 (or equivalently by their renormalized density 
matrix 7ren — 1 — -f/2). The energy of any such Hartree-Fock state reads 

(H^')=£^f(7-//2) + C 

where the constant C diverges for infinite volume, and 

(4) £:^F(7rcn)=VFtr(i?°7,en)+ / V^(x)p^_ (x) dx 

1 [[ P....i^)P...Mdxdy^^-!! "^^^^^P^^dxdy. 



2 77r2xk2 \x-y\ 2 77r2xm2 \x-y\ 

The reader can recognize the well-known Hartree-Fock energy |38| . but applied 
to the renormalized density matrix 7ron = 7 ^ ^/2 instead of the usual density 
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matrix 7. This is of course a consequence of our choice of a charge-conjugation 
invariant density operator p{x) in ^. The last two terms of the first hue of @ 
are, respectively, the kinetic energy and the interaction energy of the electrons with 
the external potential V . In the second line appear, respectively, the so-called direct 
and exchange terms. In Relativistic Density Functional Theory Tli, the latter is 
approximated by a function of p-^ and its derivatives only, a procedure which we 
shall not follow here. 

One of the main goals of this paper is to construct ground states for the Hartree- 
Fock energy Q describing electrons in graphene. This energy is not bounded from 
below and it is not well defined as such, because these states always have infinitely 
many electrons. But following [331 123 122] we shall see in the next section that it is 
possible to construct ground states by using a thermodynamic limit procedure. 

In this section we have considered generalized (mixed) Hartrce-Fock states, 
whose density matrix 7 only fulfills the condition < 7 < ^- This technique 
was first proposed by Lieb |34| and it is very convenient when proving existence 
results for ground states. In view of a variational principle from [34) . Hartree-Fock 
ground states are always pure in the presence of repulsive interactions, i.e. their 
density matrix is automatically a projection in the end. 



3. Fermi velocity enhancement in Hartree-Fock graphene 

In this section we consider a graphene sheet without any external field, 1/ = 0, 
and we investigate the effect of the Coulomb interactions among the electrons. 
In mean-field theory it is well-known that the effective Fermi dispersion relation 
becomes singular at 0. This enhancement of the Fermi velocity has already been 
remarked in [201 IS21 ESI S] ■ Our main contribution in this section is the rigorous 
proof that the so-obtained state is actually the true ground state of the system. 
The method thereby follows that of [29) . 

If the electrostatic interactions between particles are neglected and the system is 
confined to a box with periodic boundary conditions, it is obvious that the unique 
minimizer is the (non-interacting) free Dirac sea, which converges in the thermo- 
dynamic limit to the infinite-volume (non-interacting) free Dirac sea. The latter is 
an infinite Hartree-Fock state containing all the negative energy electrons (in ac- 
cordance with the old Dirac picture [319]), whose density matrix and renormalized 
density matrix are, respectively, given by 

(5) P"=1(I?"<0) and 70^„ = -^=fO-^. 

When interactions are taken into account, the free Dirac sea changes but it stills 
remains translation invariant. The latter was rigorously proved in the (massive) 
3D case in [55], for a < 4/7r. The same result will be true here. 

The energy per unit volume of a translation- invariant state 7ren = fion{p) is 
given by 



— ^ i;f / Trc2 (cr ■ p /rcn(p)) -7: 

^1^) \ Jb(o,A} 2 7r2 



(6) J'(7i.cn) = 7^^ ( VF I Trc2((T •p/i.cn(p)) - ^ / -^^^ dx 



This energy is bounded from below provided an ultraviolet cut-off A is inserted. In 
the context of graphene, the ultraviolet cut-off A mimics the presence of the carbon 
lattice in graphene. Its physical value is A ~ 0.1 A , see [48]. The state 7i.on is a 
multiplication operator by the 2x2 matrix /rcn(p) in Fourier space, supported in 
the ball of radius A. Its kernel in x space is then given by (— 27r)^^/i.cn(a:: — y) where 
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/ren IS the Fouricr inverse of /rcn- The density of charge of any such translation- 
invariant state is found to be constant in space: 

p^_ = (27r)-i trc2 Ae„(0) = (27r)-2 f trc^ (/rc„(p)) dp. 

Jb{oa) 

The second term in the energy ^ is the exchange term per unit volume. We have 
assumed that our translation-invariant state 7ron has no density of charge and thus 
there is no direct term in the energy. We will verify below that, indeed, p-^^^^ = for 
the minimizer. On physical ground it is clear why this must hold since the Coulomb 
energy of a constant density of charge is not proportional to the volume unless it 
vanishes identically. Recalling that 71.cn — 7 — -^/2, the constraint < 7 < / then 
takes the form 

(7) - ^ < Unip) < ^ for a.e. \p\ < A, 

where 1^2 is the 2x2 identity matrix. 

Remark 3.1. Let us remark that adding the ultraviolet cut-off A is equivalent to 
replacing the one-particle Hilbert space L^(R^,C^) by the Hilbert space 

(8) SJA:={ipeL^{m.^,C^) : supp(<^) C B(0, A)}. 

In this section we will show that the non-interacting state 7°gjj, defined in ([5]) 
and which consists in filling all the negative energies of the free Dirac operator, 
is the unique ground state of the interacting energy per unit volume ©. This 
surprising fact only occurs because of the absence of a mass. It is not true for 
massive particles, for which the interacting ground state depends in a nonlinear 
manner on interactions [37l [29] . 

The renormalized density matrix j^^^ is the multiplication operator in the Fourier 
domain by the matrix 

2bl 

Because the Pauli matrices are trace-less, the charge density of this state vanishes, 
i.e., PjO = 0, as was announced before. The mean-field (Fock) operator of this 
state is given by 

(9) v'>=vfD°^^^^^-^. 
^^ 27r|x-y| 

It is nothing else but the derivative of the energy (|6]) at /"^n • The second term on 
the right is the exchange term. The following gives the formula of 2?" in Fourier 
space. 

Lemma 3.2 (Effective velocity of graphene). With /°en(p) — """ ' P/(2|p|); the 
mean-field translation-invariant operator ^ can be written as 

(10) ^°(p) == t^cff (p) cr ■ P, where v^sip) := wf + 5 rr 



/r"cn(p) = -Tmr = l(-oo,0)('^ ' p) 



\P\ 



and 

(11) 



giR)^^rr ^"^^ rdrdO. 

2n Jo Jo Vr^-2rcose + l 



The function g is increasing on [l,oo). It satisfies 

on _ 1 
3(1) = ^^ ^ 0.1324 
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where G ~ X]n>o(~-'^)"(^" '^ 1)^^ — 0.9160 is Catalan's constant and 

We see that the effective Fermi velocity Ves{p) is logarithmically divergent at 
p = 0, i.e. 

Vcs{p) = - log — + O(l)|p|^o- 

This is the well-known velocity enhancement mentioned in the title of the section 
(for comparison see, e.g., [?1, Eq. (220)]). Here the 0(1) is independent of A. 

Remark 3.3. The logarithmic divergence is sometimes called the Kohn anomaly. 
It has the effect of reducing the density of states near the Dirac energy [IT] . 

Using that g{A/\p\) > g{l), we see that 

\V'ip)\>{v^ + gil)W{p)\, 

an inequality that will play an important role later when we will show that j^^^ 
is the unique minimizer in the absence of external potentials. But before we shall 
state the proof of Lemma [ 



Proof. Using that the Fourier transform of |a;| ^ is exactly \k\ ^ in 2D, we can 
write the translation-invariant operator 2?" defined in ^ in Fourier space as 



(12) WFCr-p-P"(p) = — / '-^S^dk^- — \p\cT- ,, , dfc 



with ujk '.= k/\k\. It is clear that the vector 

Wfc 



/ 



|fc|<A/|p| \k- LOp 



■dk 



is co-linear to p. Hence we can also write 



dk = LOr, / TT-^^ r dk 



'|fe|<A/|p| \k-ojp\ "'|fe|<A/|p| \k-ojp\ 

which leads to V^i^p) — Vctf{p)cr ■ p with Wcff (p) as in (fTO|) and 

9iR) = 7- / TT—1 dk 



47r./|fe|<fl |fc-Wp, 

-rdrdO 



47r 7o 7o ^/r'^ ~2rcos9 + 1 

I /■'r/2 /■« / 1 1 



, , cos9 , -. , rdrdO. 

27r Jo Jo \ V r-2 - 2r cos 6* + 1 Vr^ + 2r cos 6* + 1 / 

Since the integrand is non- negative, it is now clear that g is increasing on [1,cxd). 
For large R we have 

giR) ^ '^r\osHdO='^^. 

^^ ' R-^oo TT Jo 4 

For the value of (7(1), we integrate first in r and obtain 



/ 



\/r'^ — 2rcos0 -f 1 



dr = -l + cos(6i) log (l + sin"i (61/2)) +2sin(6i/2) 
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and we therefore get 

g(l) = -!- / (cos2(6l)log(l + sin~\6l/2)) +2cos(6')sin(6'/2))d6l 

= _L /" cos^{e)\og{l + sm{e/2))de-^ f cos^{9)\og{sm{9/2)) d9~^. 
zTT Jo 27r Jq Stt 

The result follows by explicit integration, using that 

^7r/4 n'''/4: 

= 2/ log(2 cos(6')) de = - log(2 sin(6l)) dO. 

Jo Jo 

This concludes the proof of the lemma. D 

Because V^ is equal to the original D° multiplied by vp + g{A/\p\) > g{l) > 0, 
we have 

7r"c„ = ~^j^ = l(-oo,0)P°) - 2" 

In other words, the non-interacting free Dirac sea solves the nonlinear equation of 
the interacting system. This is in stark contrast with the results of [37l [29] in which 
the interacting Dirac sea was found to be very different from the non-interacting 
one, as we have already mentioned. With Lemma 13.21 at hand, we are now able to 
prove that 7j'(,„ is indeed the global minimizer of T. 

Theorem 1 (Ground state of free graphene). Fix any ultraviolet cut-off A > 0. // 
(13) vj.>-^yj^~gil)^ 2.0560, 

then the minimization problem 

min I J"(7,.en) : 7rcn = /rcn(p), -li3(0,A)/2 < /rcn < 1b(0,A)/2: P7rc„=0| 

with J- defined in ©, admits the unique minimizer 

We note that the case of graphene in which vp ~ 2.2 is covered. 

Proof. The proof is exactly the same as in the massive case |29j , since the argument 
of [29] relies on the fact that the Coulomb potential can be estimated by the effective 
mean-field operator, which does not require a positive mass. 
Going back to Formula ©, we write 

-^(7) - -^(7?en) -jJ^i [ Trc. (V^ip) ifip) - fLm 

1 f \if - fLJi^)\' 



, , dx 

2 

Using that -l(cr • p < 0) < f{p) - /°c„(p) < t{a ■ p > 0) exactly as in [29j, we see 
that 

Trc.(l?"(p) ifip) - /°n(p))) > Trc.(|P°(p)| (/(p) - /°„(p))') 

> (vp + ff(l))Trc2 {\p\ ifip) - Aen(p))') 

for a.e. p E 5(0, A). Hence we conclude that 



Hi)-HiL)>77^{ivF+g{i)) f |(-A)V4f|2_i /■ \I\dx 

7r2 ' 2 JTg^l \x\ 



1 
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with F := f — /^'Jjjj. Kato's inequality in 2D tells us that 

(») A<i['"^''-^ 



|a;| - 2r(3/4)2^ 

where the constant is optimal (311 37] (see also [361 Lemma 8.2]). We deduce that 
^il) > -^(Trcn); providcd vf satisfies the inequality P^ . Since there is no other 
optimizer than for Kato's inequality, even the equality in (J13p is covered. D 

In summary, we have proved that 7^g„ — P^ — 1/2 is the unique minimizer of the 
energy per unit volume when V = 0. Arguing exactly as in [53], it is then possible 
to prove that j^^^ is also the thermodynamic limit of the true ground states of the 
Hartree-Fock energy, without the translation-invariance ansatz. The proof is even 
much easier than in [35] since 7°g„ is known exactly and solves the self-consistent 
equation in a box as well. Instead of pursuing this route in detail, we accept that 
7°gjj is the actual free Dirac sea, and we now study local perturbations of it, in the 
spiritof [S1I151I2S]- 

4. Ground states of Hartree-Fock graphene in local external 

potentials 

Let us now come back to the Hartree-Fock energy (|3|) and assume that the 
external field V does not vanish. We will always make the assumption that V is 
local in a sense to be made precise below, which puts us in a situation where we 
can think of the sought-after Hartree-Fock ground state as a local perturbation of 
the free Dirac sea. 

We consider any Hartree-Fock state described by its density matrix 7 (or equiv- 
alently by its renormalized density matrix 7rcn — 1 ~ -^/2) and which we as- 
sume to be "sufficiently close" to P^. The infinite volume Hartree-Fock energy 
of 7rcn is of course infinite, it is proportional to the volume exactly like the one of 
7ren — P'l ~ ^/2. Howcvcr we can, at least formally, subtract the (infinite) con- 
stant iSHF(7rcn) — '^HF(7rcn) ^^'^ obtain a perfectly well-defined energy. A formal 
computation yields 

^HF(7rcn) — ^HF(7rcn) — ^BDf{Q)j 

where 

Q = 7ron - 7"cn = 1 - P- 

and where S^^p is the so-called Bogoliubov-Dirac-Fock energy, formally defined by 

(15) f^DF(g) = tr(P"Q) + / V{x)pq{x) dx + l ff ^^MM^lrf^ dy 

Jk2 2JJr2xk2 \x-y\ 

' ff ^Q^^^y^dxdy. 



2JVxE2 \x-y\ 

Again the energy functional looks like the usual Hartree-Fock energy, with the 
difference that !)'•' now appears instead of D^ and that it is applied to the operator 
Q which is a difference of two Hartree-Fock density matrices. The operator Q 
satisfies the constraint 

-P- <0 < 1-P° =P+. 

Remark 4.1. To our knowledge, the idea of subtracting the infinite energy of the 
free Dirac sea in order to get a bounded below energy, was used for the first time 
in [29,. This was generalized to positive temperature in [24) . In previous works [5] 
[53] dealing with the Hartree-Fock approximation of QED, another justification 
based on normal ordering was employed. 
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Remark 4.2. Let us mention that the perturbed state 7 can always be seen as a 
Bogohubov rotation of the free Dirac sea in its Fock representation, which is why 
Chaix and Iracane used the name 'Bogohubov' for the energy (TTSl) . We could as well 
call it a relative Hartree-Fock energy but we prefer to keep the name Bogoliubov- 
Dirac-Fock (BDF) for historical reasons. 

Our tasks in this section are then to prove that: 

(a) Q = is the unique minimizer of fgup for V = 0, which is a "local" 
version of the fact that the free Dirac sea P^l is the unique ground state of 
the system without external field; 

(b) if y 7^ 0, then there exists a ground state for fg^p which solves the self- 
consistent Hartree-Fock equation. 

Before turning to these problems, we however have to properly define the BDF 
energy. It is now well understood that the Hartree-Fock ground state 7 of an 
infinite Coulomb system can in general behave badly. Usually Q = 7 — P- is not 
a trace-class operator and its density is sometimes not in L^, see [22j [5]. It has 
long-range oscillations which are not integrable at infinity in some cases. For these 
reasons, it is not fully obvious to give a clear meaning to the BDF energy P3|) 
and to find a suitable class of states in which minimizers will be. Following ideas 
from |25) . we introduce the correct functional analysis setting in the next section. 

4.1. Function spaces and definition of the density. Given an operator Q, we 
define Q'^'^ := P^QP^, where e, e' £ {±}. Our starting point is the remark that, for 
a nice-enough operator Q (say finite rank), 

tr(pOQ) = tr (pO|(Q++ - Q—) > tr \V°\Q\ 

Here we have used that P° commutes with P'^ and that 

~P^<Q<Pl ^^ Q^<Q++-Q—, 

as was remarked first in [T]. We see that a state will have a finite relative kinetic 
energy when |I?°|-'^/^Q^^|2?°r^^ a-re trace-class, but we cannot gain any other in- 
formation on Q^^ than jP^j^/^Q^^ being Hilbert-Schmidt. Thus we shall assume 
that 

(16) |I?°|l/2Q±±|2?0|l/2g@l ^^^ g±T|pO|l/2^g2 

where 6^ denotes the usual p-th Schatten space (6^ and (3^ are respectively the 
spaces of trace-class and Hilbert-Schmidt operators). This enables us to properly 
define 

tT{V°Q) := tr (|P°|i/2(Q++ - Q—)\V'>\'/^) . 

Our next task is to give a clear definition of the density pg under the very weak 
assumptions (J16p on Q. To simplify our exposition we introduce the Banach space 

where we recall from ([5]) that Joa is the Hilbert space of L^ functions with compact 
support in B{0,A) in the Fourier domain. The set of bounded operators on this 
space is then denoted by B{^\). We also introduce 

i- := {q e S(i3A) : Q* = Q, |P°|1/2q±±|2?0|1/2 ^ @i^ g|pO| ^ g2| ^ -^^ 

since jP"! > {vp + g{l))\D^\ — {vp + g{l))\p\, as we have shown before. For what 
follows it will be convenient to work with states in X. By doing so, we ignore the 
logarithmic divergence at of 'D°(p), which is an additional information for us to 
be used in due time. Let us remark that X contains elements which are not even 
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compact. This makes the mathematics more involved than in the situation where 
the particles have nonzero mass. 

Next, we shall show that any Q G X indeed has a well defined density pg. We 
start by remarking that pg is locally well defined. 

Lemma 4.3 (Definition of the density pq). Any Q (z X is locally trace class, i.e., 
for all x{x) G L°°{M.'^) with compact support, xQx G 6^ ■ Moreover, the density pQ 
is inL°°(R2). 

Proof. Due to the cut-off we have 

xQx = xHaQHax 
where Ha — l(|p| < A). Since x € L^, then clearly xIIa G S^, such that together 
with the boundedness of Q we get xQx S ®^i ^s stated. This proves that pg is a 
well defined function in L^^^. To see that it is actually uniformly bounded, we use 
that Q is self-adjoint and the ultraviolet cut-off A, to infer — HQHIIa < Q < HQHIIa. 
On the diagonal we get \pq{x)\ < WQWpn^ix) = ||(9||AV(47r). D 

Remark 4.4. Any bounded sequence (Qn) in X has a weakly—* convergent sub- 
sequence, Qn^. — ^ Q in the sense that 

tr(AQ„J ^ tr(AQ), VA e 6\ 

tr(i^|p|i/2Q±±b|i/2) ^ tr(i^b|i/2g±±|p|i/2)^ VK compact, 

and 

tr(BQ„, |p|i/2) ^ tr(BQ|p|i/2), VS e S^. 

Using that xIIa G 6^ when x & L"^ , it is then elementary to verify that xQnkX ~^ 
xQx strongly in the trace-class. This implies that pq^ -^ pq strongly in L^^^, 
hence strongly in Lf^^CM."^) for all 1 < p < oo, by interpolation. 

Next, we introduce the so-called Coulomb space: 

C:=\ip : D(ip, if) := 27r / MML dk < oo 
{ Jw^ \k\ 

which is the natural energy space for the density pq. If we decompose the density 

of Q into 

PQ ^ PQ++ + Pq-- + Pq+- + Pq-+^ 

then we can deduce the following properties for the elements of X. 

Lemma 4.5 (The density is in C). Assume that Q E X is such that — P'^ < Q < ^+- 



les 



Then pq±± £ LP{M. ) for 3/2 < p < oo and pq±T & C, which particularly impli 
that Pq £ C. 

Proof. Since < ±Q^^ < 1, the Lieb-Thirring inequality in 2D for the relativistic 
kinetic energy [36] immediately yields 

±tr(b|Q±±)>c/ \pq±±{x)f^^dx. 

Since we already know that pq±± € L°°{M.'^), pq±± G C now follows from the 
Hardy-Littlewood-Sobolev inequality in 2D, see, e.g., [35J. The second part of the 
statement follows by a duality argument, provided that we can show 

|tr^Q+-| = |(e,PQ+-)|<c(^^Jfc|||(fc)pdfc) , 
for aU ^ e C~(M2). To this end we first estimate 



treQ+-=tr(^^e^HV4QHV4)< 



pO pO 



llblQIIe- 

62 
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Using 



Trc.P°(rtP^(g) 



bikl 



we can compute 



pO pO 

11/4 "S I 



|p|l/4-|p|l/4 



e^ 



<^ ff ,.n.,2 \i + k/m-k/2\ {i + k/2) . {i-k/2) 

,2,,, f \l + LOk/2\\£~iOk/2\-{i + ujk/2)-{i-LJk/2) 



: Je2 \£ + uJk/2\-^'-^\£-uJk/2\'^'^ 

<c / i£(fc)nfc|dfc. 



Here ujk '■= k/\k\ and we have used the fact that 

1 f\e + LJk/2\ \i - Uk/2\ -{£ + Uk/2) ■ {£ - LUk/2) ,^ 

at = I 



is a finite integral, independent of the direction ujk E S^. D 

4.2. Stability of the free Dirac sea. We have shown that the density pQ is in 
the Coulomb space C whenever Q G X and — P^ < Q < P+. Thus we see that the 
direct term is well defined for all such Q with finite relative kinetic energy. We can 
now define the (free) BDF energy as 

£"bbf{Q) ■■= tr(I?°g) + Id{pq,pq) -\ff MMl^dxdy. 

where we recall that 



D{p,p')-= -. T-dxdy^2T: — dk 

JJm^xr^ \x~y\ Jr2 \k\ 

is the so-called Coulomb scalar product. 

The following lemma shows that the exchange term is also well defined and that 
the BDF energy fg^p is non-negative with Q = being its unique minimizer. 
Recalling that fgop is the relative energy counted with respect to the free state 
P^ and that Q = 7 — P° , this consequently shows that the free ground state P° 
is stable under local deformations. Here 'local' refers to perturbations such that 
Q € (Y but not necessarily small in norm. 

Lemma 4.6 (Stability of P^). For any fixed ultraviolet cut-off A, the mapping 
Q H^ ^bdf(Q) *'5 ™6^^ defined and continuous on X . If vp satisfies 

ir(i/4)2 

then we have 

for all -pO <Q <Pl with |X'"|i/2g±±|X)0|i/2 ^ gi^ Furthermore, E^j^piQ) = 
if and only if Q = 0. 
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Proof. Using Kato's inequality (|14l) we infer, following [1 and similarly as in the 
proof of Theorem [I] 

2 yy \j:- y\ ^^^2^ ^ 4 r(3/4)2 ^'^^P^^ > " 4r(3/4)2(z;F + .g(l)) *'^l^ l'^ ^' 

which is well-defined due to our assumption that Qjl?''!^/^ S 6^. In addition, since 
-P- <Q <Pl, we have trd^^lQ^) < tT{V"Q) and therefore 

(18) i [[^^^d.dy< ^^^,l^:/'^' , tr(pOQ). 



2 77 |2:-y| ^-4r(3/4)2(vF+5(l)) 

In other words, if |X'0|1/2q±±|X)0|1/2 ^ @i and if the Fermi velocity wp satisfies (|T7)) . 
the exchange term is controlled by the kinetic energy. We have already seen that 
PQ G C. Using then that D{pq,pq) > 0, we see that S^DpiQ) — ^ ^^^ ^^^^ Q = 
is the unique minimizer of ^bdf' ^^ stated. D 

4.3. Existence of minimizers for graphene in an external field. In the next 
step, we shall submit our graphene sheet to an external electrostatic field of the 
form 

with 1/ the density of charge of the defect. The corresponding BDF energy now 
reads 

Using our estimate (TTSl) on the exchange term and that D{-,-) defines a scalar 
product, we immediately get the lower bound 

fV(Q)>"^^('^,^), 

provided that vp satisfies (fT7|) . Therefore the BDF energy is bounded from below 
when ly £ C which is our way of measuring the locality of the potential V . This 
enables us to consider the minimization problem 

(19) E^j,p -.^ M [s^^^iQ) : 

Q e S(Joa), -p°<Q^Q*< Pi, |p°|i/2g±±|i?0|i/2 g gij. 

We shall prove the existence of a corresponding minimizer, which is the Hartree- 
Fock ground state of interacting graphene in the presence of defects. This existence 
result is non-trivial for the simple reason that we have no mass. 

Theorem 2 (Existence of a ground state for infinite volume graphene with defects). 
Fix A > and let v-p be such that 

ir(l/4)2 

For any v in the Coulomb space C, the problem (|19p admits at least one minimizer 
Q, satisfying the self- consistent equation 

(21) Q + pO=i^('pO_^ 1 ^^^,^_Q(^ 



Fl Fl \x~y\ 

where I — (— oo,0) or I ^ (— oo,0]. Equivalently, with J — Q + P^ denoting the 
density matrix of the optimal HF state, we have 

^ f , ^, 1 1 h-I/2)(x,y)\ 

7-lx(.F^-(-^V)-.*^+p,_./.*^-l^-^^A^j. 
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Remark that HF ground states of graphene in the presence of a defect can be 
chosen pure and with the Fermi level either filled or unfilled completely, as is usual 
for Hartrce-Fock theories. This follows from Lieb's variational principle f3l[ and 
the no-unfilled shell theorem of Bach, Lieb, Loss and Solovej [2j. Note the strict 
inequality in (I^Ul) . There are probably also ground states in the case of equality, 
but we shall not consider this case for the sake of simplicity. 

The proof of theorem[5]is a bit long and it will be given in Section [S] below. The 
method is similar to the one in [55], but several modifications are needed due to 
the absence of the mass. One important additional input is a localization estimate 
inspired by [33] and which is detailed in the Appendix. In contrast to 25] we shall 
use that the kinetic energy has an infinite velocity for p ~ 0, which induces a better 
control of the exchange term. Our method docs not seem to apply otherwise. 

Remark 4.7. We do not know if the optimal state Q has a finite (relative) num- 
ber of particles, that is, Q^^ might be not trace-class and the charge tr((5) := 
tr((5++ -|- Q ) could be infinite. The operator Q could even be non compact in 
general, because of the absence of a gap. This would mean that the corresponding 
minimizers Q live in a Fock representation which is not equivalent to that of P^ , 
even if the relative energy is itself finite, by the Shale- Stinespring theorem [JS]. 
In |27j another minimization problem consisting in fixing the relative charge tr((5) 
was considered. Because there is no gap and tr((5) can be infinite, an analogous 
approach does not seem to make sense in our context. Bound states with finitely 
many electrons have been constructed in a projected Dirac-Fock-type model in |10) 
but there a magnetic field is used to confine the particles and create a gap. 

5. Linear response to an applied external field 
In this section we consider a small external field 

Vx = -\i^ * —r, A < 1 

m 

and discuss the linear response of graphene within our Hartree-Fock theory. To this 
end, we denote by Q\ a chosen minimizer for each A and note that £^^p{Q\) < 0, 
which is seen by using Q = as a trial state. From this we deduce that 



2 77r2xr2 \x-y\ 2 ^2 



II 2 



c ■ 



Assuming the strict inequality 

1 £(1/4)2 

4r(3/4)2 
which allows to control the exchange term, in view of (|18p . we deduce that 



(22) -^>774t77TI-5(1), 



tT\V"\Ql+ ff 



\Qxix,y)\' 
\x-y\ 



dx dy + II Pa — '^'^llc = 0{X^ 



with p\ := pq^. This confirms that Q\ is of order A in X. 

Recall from Theorem [2] that Qx satisfies the self-consistent equation 

Qx - 1(-oo,0)(2?qJ - P" = 1(_oo,o)(2?qJ - l(_oo,0)(2?"). 

For simplicity we assume here that ker(2?Q^) — {0} for all A. Denoting 
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with (pA = Ai' * |iy - Pa * ||y, and Rx = Q\{x,y)/\x - y\, we can write, using a 
formula of Kato, that 



^^-i 



1 1 



Vq^ +iri X>o + iri 
In view of the resolvent formula we deduce that 



drj. 



Q\ = Ql.D + Qi,x — TT' 

Ztt 



1 . N 1 

-{(fix + Rx)- 



dr] + 0{\^ 



The remainder O(A^) has to be estimated carefully, but for convenience we remain 
formal in this discussion. Because ipx and Rx are themselves affine in Qx, the first 
order term for Qx is obtained by inverting these linear maps. There is no simple 
expression for it. However, following [301 US], it is possible to compute explicitly 
the density coming from the term involving (px, as will be explained now: 

We look at the density associated with the operator stemming from the direct 
term 

2n J^ V" + irj V" + irj 
In Fourier variables this reads 

\ [^ I 1 

Using 

1 /■°° 1 1 1 

M{p,q) = - T,a/ \ , ■ T,of \ , ■ '^^= T?f \ , j^t j '7--a;pcr-cjg-l), 

TT Joo ^ {P) + iV ^ {q) + ^V E[p)+E{q) 

with 

E{p)^\p\(v^+g( — 

we can write 

Qi,D{p,q) = —0{p-q)M(j>,q). 

Using 'D^{p)/\E{p)\ = a ■ ujp, we see that the corresponding density reads 
Pi.oik) ^ TT Trc2(3i,_D(p + k/2,p~ k/2)dp 

^^ J \p+k/2\<A,\p-k/2\<A 

= -l^A(fc)|fc|i?(fc) = {Oik) - pxik))Bik), 

ZTT 

where we have used that in 2D 

m \k\ 

and we denote 

„,,^ l_ f ip + k/2)-ip^k/2)-\p + k/2\\p-k/2\ 

^ > \k\27r 7!p+^/2||A \p + k/2\ \p - k/2\{{E{p + k/2) + E{p - k/2)) ^ 

\p—k/2\<A 

1 f ^{p + uJk/2) ■ {p - uJk/2) + \p + ^fc/2| | p - a;fc/2| 

X 



27r yiP+-./2|<A/|fc| \p + tjfc/2| |p - ujkl2\ 

|p-i.jfc/2|<A/|fe| 

1 

v 

b + c.,/2|(t;F+g( |,||^^^^^,| ))+|p-^fc/2|(^F+g( |,||,l,/,| )) 



dp. 



Remark 5.1. In the case of three spatial dimensions there is a similar function 
B{k) playing an important role. In this case, the value S(0) > is logarithmically 
divergent with respect to A, which is the reason for the requirement of charge 
rcnormalization, giving rise to the Uehling potential, see [501 [^ [^ I23j . 
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In Fourier space, we can write the self-consistent equation as in |26| 
pl(fc) = B{k){u{k) - Mk)) + Pi,x{k) + O(A'), 
where pi^x — 0{X) is the first-order density coming from the exchange term 

1 f°° I 1 

Qux ^ TT / ^0 I ■ -^^ -no I ■ ^^^ 
ZTT J^ T>^ 4- irj V" + irj 



We see that 



«<')^^Tlil^<« + ^li-«'^'^ 



It is therefore natural to ask about the behavior of B{k) for low momenta. If we 
neglect the density pi,x stemming from the exchange term, this will determine the 
decay in x space of the Coulomb potential of the (first order) polarized graphene 
in presence of the external density v. 

In order to answer this question, we simplify the expression of the function B{k). 
First we remark that B(k) is obviously radial, hence we can take u)k = £i '■= e, 
such that 

(23) Bik)- ' f b + e/2|b-e/2|-p^ + l/4. 



2Tr J\p+e/2\<A/\k\ \p + e/2\\p-e/2\ 

\p~e/2\<A/\k\ 



dp. 



\p + e/2| [vF + 9 { \k\\X/2\ ) } + \P- e/2| [vF + 9 (^b^) } 
As in [531 [35], we use the following change of variables, 

_ \p + e/2\~\p~e/2\ _ \p + e/2\ + \p - e/2\ 

'" 2 ' ^ 2 ' 

Denoting p — [x, y), this reads 



^{x + l/2)^ + y^~^ix~l/2y+y^ 
V =— 



^{x + l/2y + y^ + ^{x~l/2y+y^ 



The corresponding Jacobian is 
d{v,w) 



d{x,y) 



y \P2 



2 V(x -I- 1/2)2 + y2^(3, _ 1/2)2 + y2 2\p + e/2\ \p - e/2| ' 
We collect the following relations 

2w^ + 2v^ = {w + w)2 + {w- v)^ = 2(/ + 1/4), 

Avw = {w + vY — (ui — v)"^ — 2p ■ e = 2pi = 2x, 
and 



\y 



= |P2 1 = Vp^ - (P ■ e)^ = V^^ - 1/4 - 4^2 (w2 - 1/4) 



= 2v/w2 - 1/4^/1/4 -w2. 

Observe that 

^^|(p + e/2) + (p-e/2)|^ 

2 ' 

and 

|.| < \(P + ^/^)-(P--/^)\ ^ 1/2. 

Our constraints onp, i.e., \p+ujk/2\ < A/|fc| and \p—uJk/2\ < A/\k\ can be expressed 
in terms oi w,v as follows 

w < A/\k\ and \v\ < w - A/\k\, 
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such that 

1/2 < W7 < A/|fc|, 

\v\ < niin{l/2,A/|fc|-i(;}. 
Changing variables wc find 

|p + e/2||p-e/2|-p2 + i/4 



dx dy 

-,„2 _ 1 M^ _|_ 1 M 9lT)4-p/9lln_ p/9l 

dv dw 



\p + e/2\\p-e/2\ 

_ (w + w;)(u;- w)- (i;2 +u)2 - l/4) + l/4 2\p + e/2\\p - e/2\ 



\p + e/2\\p-e/2\ 



'^'"-^'^ ' ' ^y/y^^.^--.. 



^w2- 1/4^1/4 -w2 ^w2 - 1/4 

On the other hand we can express 

\p + ujk/2\ \vF+g( TTTp- -7 ) i + b - i^k/2\ \vF+g 



\k\\p + Uk/2\Jl ' '" -/ I ^ - ■ --"WkWp-ujkm 
{v + w) ivp + g [ —— I > + {w - v) <vf + g 



\k\{v + w)J j \ ■ \|fc|(w — u) 

to arrive at 



A/|fc| /.min{l/2,A/|fe|-u.} /l /A_„2 



2 /-A/lfcl |.min{l/2,A/|A;|-u.} ^1/4 



(24) B(fc) = - dw dv 



v^ 



TT 7i/2 7o ^u;2 - 1/4 

1 



:X 



S^Fti; + (« + ti')g (^7^) + (^i' - v)g (^^) 

Remark 5.2 (The no-exchange case). If we discard the exchange term, then we 
have exactly the same calculation with g replaced by everywhere. In this case the 
linear response involves the modified function 

1 /-A/I*:! pmin{l/2,A/\k\-w} /l /4 _ „2 

(25) i?0(fc) = — / dw ^ )'^ dv. 

TI'WF Jl/2 Jo Wy'w'^ - 1/4 

Let now w — t/2, and v = {cos9)/2, such that < cos0 < min{l,2A/|fc| — t}. 
Then 

1 f2A/\k\ (.arccos[min{l,2A/|fc|-t}] -, 

B°(k) = / dt <.\n^nM 

Since 2A/|fc| — t < 1 is equivalent to i > 2A/|fc| — 1 we can decompose it into two 
integrals 

1 / /•2A/|fe|-l ^7r/2 , 

B°(k) = / dt ^\n^flM 

27rvF\Ji Jo t^/W^l 

2A/|fe| .x/2 ^ \ 

dt / —j^=^ sin^ QdQ 

2A/|*;|-1 Jarccos[2A/|fe|-t] t\J t"^ - 1 j 

--arccos(l/(2A/|fc|-l))H / dt- 



2'KVY 4 2'KVY J2KI\k\-\ 4:ty/t^ - 1 



X ( —2 arccos ( 2-r—r — t ] + n + sin 
|fc| 



—4r' 



This immediately shows that 



lim B°{k) ^ 



k^o 16wf ' 

which obviously depends on v-p in contrast to what we will find for B{k) below. 
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Let US now come back to the function B{k). The following lemma says that the 
function B{k) vanishes at A; = 0. However, it only vanishes logarithmically. 

Lemma 5.3. We have that 

B{k) 



" 4 log [{t 
Proof. In ([M)) . we split the integral in w as follows: 

A/|fc| .log2(A/|fe|) .A/tfct 

dw — I dw + / dw. 

1/2 Jl/2 Jlog2(A/|fc|) 

By doing so we can write B{k) = Bi{k) + -82 (fc) with an obvious definition. For 
w e [log (A/j/cl), A/|fc|] we just use that g > g{l) and we obtain 

1 /-A/lfel ^min{l/2,A/|fc|-«;} /i 74 _ ,,2 

B2{k)<— -TTT / dw dv^ ' 



"■(^'F +5(1)) Jlog2(A/|fc|) Jo w^w"^ - 1/4 

Using the computation of Remark l5.2l we see that the right hand side behaves like 
1 /tt / 1 W 1 



8(^F+5(1)) V2 Vlog'(A/|fc|);; IfcKo 8(t;F + <?(!)) log2(A/|fc|)" 

In particular, 

limlog(A/|fc|)B2(A:) = 0. 

fc— J-O 

Now for w G [1/2, log (A/|fc|)], we can safely expand the terms involving g. Indeed, 
we have w ±v < log (A/|/i;|) + 1/2, hence 

A A 

> 



\k\{w±v) '- |fc|(log2(A/|fc|) + l/2) fe^o 
By the dominated convergence theorem we find that 



A ^log2(A/|fc|) .mm{l/2,A/|fe|-u,} , /TTJ" 

limlog(A/|fc|)Bi(fc) == lim - / dw dv ■ ^ ' 



fc^O fe^O TT Ji/2 Jo Wyjw'^ - 1/4 

The right hand side is again similar to B^ (fc) computed in Remark 15.21 and it is 
equal to 7r/4. D 

Note that, in contrast to the no-exchange function B^ which has the wp-dependent 
finite limit 7r/(16fF) at fc = 0, the true function B{k) tends to zero and its behavior 
on first order is universal, it does not depend on v-p. As we have explained, if we 
neglect the exchange density coming from R\^ then we find that in first order in A 

p^(k) ~ A ^^^],M k) ~ A l^^v{k). 

^'^^ ' A«i 1 + B{k) ^ ' A«i 4 ,^p, {a} 

^ ' \k\«\ * iog ifci , 



|fc|«i ^'^^\\k\ 

That this density vanishes at fc = means that the response of the graphene sheet 
in the presence of v is essentially neutral. This is stark contrast with the massive 
3D case, in which there is always (partial) screening, i.e., -B(O) > 0. However the 
fact that here Bik) — >■ only logarithmically creates some long range oscillations 
in X space which induce some weak screening effects. Indeed, if we compute the 
first-order polarization charge of graphene in a ball of radius K (using a smooth 
localization function x), we get 

PQix)x{x/R)dx^X [ ^^^^^ ^ xiRk)R'dk ^ X-''^^''' 



R2 4iogrA\^^ fl^^ 41og(i?)' 
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It goes to zero when i? — >■ oo, but very slowly. So in a finite ball of radius R we 
might have the impression that Jg pg > is not small, i.e., that the external 
density v is partially screened by the spontaneaous polarization of graphene. 

In this discussion we have neglected the first order density pi,x coming from the 
exchange term Rx. The effect of this additional term is not clear to us. 

6. Proof of Theorem [2] 

Let (Qn) be a minimizing sequence for ([T^. Using the strict inequality (PH)). it 
is easy to see from our estimates that {Qn) is uniformly bounded in X: 

(26) tr(|P°|i/2Q2|pO|i/2) <^^ 

(27) tr (|I?0|i/2(Q++ _ Q--)|p0|i/2) < c 

Therefore, up to a subsequence, there exists an element Q Cz X , such that 

Qn ^ Q weakly-* in X, 

similarly as in Remark 14.41 This implies in particular that pg^ -^ pq weakly in C 
and strongly in L\^^ (R^ ) , by Lemmas 14.31 and 14.51 Since we also have 

{ip,Qn(p') ->■ {f,Q(p'), y(p,ip' (ES)a, 

it is clear that the constraint —P^<Q< P^ is satisfied. Therefore Q is an 
admissible state for the minimization problem (J19p . 



Step 1: Our goal is to show that the energy is lower semi-continuous, i.e. 

(28) liminf5^(Q„) >£^(Q). 

n— >-C30 

Observe that the function 

p^ D{p-v,p-v) 
is lower semi-continuous. Therefore 

liminf D(pQ„ - v.pq^ -v)> D{pQ -v,pQ-v), 

n— ^OO 

using the simple fact that pq^^ -^ pg weakly in the Coulomb norm C. Since 

tT\V°\iQ++-Q--) 

= tr {\V'\ - \p\{vp + 5(1)))(Q++ - Q--) + [VF + ff(l)) tr |p|(Q++ - Q--), 

and I'D'^l — \p\{vf + g{l)) > 0, we can use Fatou's Lemma for trace-class operators 
to obtain 

tr {\V°\ - \p\{vp + 5(1)))(Q++ - Q--) > tr [|2?°| - \p\{vf + 5(1))] (Q++ - Q^), 
and (|28p will be achieved by showing 

(29) liminf ftr b|(Q++ - Q--) | // \9lli^dxdy 

" V 2(wF+5(l)) JJr2xr2 \x-y\ 

>trb|(Q++-Q-)-— ^-^ // ^^^dxdy. 



2{vF+gil)) JJr2x^2 \x-y 
This inequality will now be proved in Step 2. 

Step 2: The proof of (E^ will heavily rely on (P7)) . i.e., the uniform boundedness 
of Qn in the energy norm 2?°. The first ingredient of the proof is to split the space 
R^ into a region close to the defect v, and another one far away, and to show that 
the energy essentially localizes. A similar idea was used in [25] but because of the 
absence of a mass, we have to argue differently. 
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Let US consider two real functions XtV ^ (^""([OiOo); [0,1]), such that x = 1 on 
[0, 1] and X = on [2, oo) and x^ + V^ = 1- We define 

XRix) = xi\x\/R), m{=^) = v{\AlR) Vx e m?. 

Now we apply the localization Lemma lA.ll given in the Appendix (and which is a 
consequence of the localization estimate obtained in [33 ) in order to obtain 

(30) tr |p|(Q++ - Q--) > tr b|xK(Q++ - 0,T")Xfl + tr b|r7fl(Q++ - Qn-)m 
-c (tr b|(g++ - Q--)f' (llVxfilli. + wyvRWh)'^' (llVxi^lli. + llV^yi^Hi,)'^' ■ 
Since 

livxfllli. = livxiii^, ||vxi?||i4 = ;^||vxiii4, 

and analogously for rji^, the uniform boundedness of tr |p|(Q^^ — Qn)i implies 

tr b|(Q++ - Q;:^) > tr \p\xR{Qt^ - Q7r)XR + tr \p\ilR{Qt^ - Qn^VR - J^, 
for an n-independent constant C. On the other hand we can write 

2 \x-y\ 

Using now 

Q++ - Q-- > Ql 

as well as Kato's inequality ([Ti)). we obtain, similar to [5^, that under the assump- 
tions on v-p 

tr mVRiQn - Qn \VR > tr IpIVrQuVR > : — ttt / / , , dxdy. 

«F+5(1) J J \x-y\ 

In order to conclude Step 2, it remains to show that 
(31) 

lim hm (trHxH(Q++-g--)x«-^^^7TTr/^^^4^^^^f^d-dy 
R^oon^cx>\ 2(vF+g{l))J \x-y\ 

-tr|p|(Q++-Q-)- im ^^d^dy, 

2{vF+g{l))J \x-y\ 

where it is important to do the limit n — > oo first and then the limit i? — )• oo. To 
this end we consider the two terms 



\x-y\ JJr^xr-^ \x-y\ 



tr|p|xi?(Q^^ -Qn )XR and / 



x|(a;)|Q„(a;,y)p 



dxdy, 



\x-y\ 

separately. First, we observe that, thanks to the cut-off 

(32) tr \p\xr{Q++ - Q--)xR = tr blxHnA(Q++ - 0,;-)nAXfl. 

As we have already seen, the operator HaXr is Hilbert-Schmidt. On the other hand 

\p\XR^A = yi ■ ([P: XrI^a + XRPT^AJ 

p I 

= -TT • (^ - *(Vx/?)nA + XflPHA 
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is Hilbert-Schmidt since p|p|^^ is bounded and (VxaJHa and xb.P^a are both 
Hilbert-Schmidt. Now we can use that Qn ^^ Q weakly-* in B{Sj\) which imphes 
in particular that iv{BQnB') — > ti{BQB') for all B, B' ^ & , and we obtain 

lim tr \p\xR{Qt^ ~ Qn')XR = tr b|Xfl(Q++ - Q'')XR- 

n— >oo 

Passing then to the limit i? — > oo gives 

lim lim tr b|xfl(Q++ - Qn~)XR = tr |p|(Q++ - Q"). 

ii— >oo ri— >oo 

It remains to prove the convergence of the exchange term in ((3T|) . Recall 

An -.^ Qn - Q ^ mX 

and that tr |I?°|Q^j is uniformly bounded. Notice first that, by the Cauchy-Schwarz 
inequality, it suffices to show 

,^ ,^ ff xii^)\M^^y)\\^,y^o 

fl^oon^oo_/_/jj,2xR2 \x-y\ 

to conclude the second step of the proof. To this end we introduce an additional 
decomposition in the infrared regime in the following form; 

An - (l(b| < e) + l(b| > e))An{ti\p\ < e) + l(b| > e)) 

^t{\p\>e)Anti\p\>e) + Bn, 

such that the first term on the right hand side 

l(b| > e)Anl{\p\ > e) := ^^ 

consists only of momenta larger than a given e > 0. Since by assumption An\p\^''^ -^ 
weakly in 6^, we see that An ^ in 6^. Due to the cut-off in Fourier space, we 
deduce that the kernel An{x, y) of An converges to zero weakly in H^{^^ x R^) for 
all s > 0, and hence strongly in L^^^{R^ x M^) for all p. Now we simply decompose 

^MEiM^.. dy = // xlix)xUy)\A^i.,y)? ^^ ^^ 
\x-y\ JJv?xv? \x-y\ 

xl{x)vlR{y)\Al{x,y)\'' ^^ ^^_ 
\x-y\ 

The first term on the right hand side converges to zero since An{x, y) — ^ strongly 
in, say, L^{B2r, x B^b) whereas x\{x)x\R{y)\x — 2/l~^ belongs to L^/^(R^ x R^). 
The second term on the right hand side is bounded by 

xli-HR{y)\A^i:^.y)\\ ^,y < c ^ ^tr|p|(^„)^ < ^. 

\x — y\ R Re Re 

It remains to consider the terms 

Bn = -B„ + S„ + _B„ 

depending on the location of the cut-off function l(|p| < e) with the only important 
point being that cither of the terms has one such IR-term l(|p| < e) on at least one 
side. Again with Kato's inequality ([H]) we bound for /c = 1, 2, 3 

xli^)m^^y)\\^,^^lf \^l^^a.dy 



y\ JJm^xR^ \x-y\ 

< 



Cmin(tr(|p|(i?^)*Bf;) , tr(|p|B^(i3^)*)) . 
Choosing on the right the term for which t{\p\ < e) hits \p\, we see that 



// 



xlix)\B^ni^,y)\ 
\x-y\ 



2 



-dxdy<Ctv{\p\ti\p\<e){Anf) 
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for fc = 1,2,3. At this point we use the infrared logarithmic divergence of the 
effective Fermi velocity to estimate 

tr(|p|l(b| < e){A^f) < ^^trP°|(A„)2 < y-^y 

log(A/e) log(A/e) 

If we take first n — >■ oo, then _R — >■ oo and finally e — >■ 0, we have concluded Step 2 
of the proof. 

To summarize, we have shown that for a minimizing sequence {Qn), 

'E'bdf = liminf £'bj3p((3„) > 6^bdf(Q) > -^bdf 

where we recall i^BDF i^ tl^*^ infimum of the BDF energy. So we deduce that 
^bdf('3) ~ ^BDF ^^^ that Q is a minimizer. 

Step 3: It remains to prove that it satisfies the self-consistent equation. This can 
be done as in [2^1 Lemma 2] with the additional problem that there is no gap in 
the spectrum of the mean-field operator. For completeness let us indicate the idea 
of the proof. Since Q is a minimizer, then for any other admissible state Q', we 
know by convexity of the constraint that 

±E''{{l-t)Q + tQ') >0, 



dt 
which implies that 

(33) tYV\Q' -Q) + D{pq-u,pq,^q) + ^ 
Let us remark that the operator 



Q{x,y){Q' - Q){x,y) 



\x-y\ 



dxdy > 



Vr 



V" + (pQ - u) * 



Q{x,y) 



\x-y\ 

is self-adjoint on the same domain as P". This follows from Rellich's theorem, 
since the two self-consistent terms are relatively bounded with respect to I?°, with 
relative bound as small as we want. For instance, we have for any if € f)\ 



iPQ -I')* TlV 



< 



LH 



1 



I ' I L4(K2 

< C||pQ-J^||c||(p||^4(R2 
<C\\PQ-V\\c\\p\'l^^ 

<C\\pQ-v\\c{e\\v\ 



M 



L^iV."^) 



L2 



2e(«F+5(l))"^" 
where we have used the Sobolev inequality ||/||i4(R2) < C |||p|^^^/||^2fR2-| and 



1 



<C 



Wp 



Li(R2) 

The argument is similar for the exchange term 

2 \ 1/2 



CD{p,p) 



1/2 



L2(R2) 



Q{x,y) 

K? \x-y\ 



j 



fiy) dy 



< 

<c 
<c 



\x-y\ 



dy 



MvW 



\x-y\ 
\x-y\ 



2 \x-y\ 

1/2 

dy] {ip,\p\ip)^^'^ 



dy 



1/2 



1/2 

dy] (ellPVII 



2e{vF+gil)) 



M 
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Taking the square and integrating with respect to x gives the infinitesimal relative 
boundedness. 
Now we choose 

(34) Q':=1(_o,^o)(2?q)-P°, 

and we claim that 



-.0/^/ ^^ , n^„ .. „ \ , ^ II Q{x,y){Q' -Q){x,y) 



tro V"{Q' -Q) + D{pQ -u,pq,^q)+^ :i^^«£i:2 ^'^ ^-' dx dy 

JjR2xR2 \X-y\ 

(35) 

= tr \Vq\{P'{Q' - Q)P' - {PY{Q' - Q){PY) 
(36) 

<~iv\DQ\{Q' ~Qf 

The first line (|55|) would just be the linearity of the trace if all the operators 
were trace-class. Because of our generalized definition of the trace, it is more 
complicated to verify ([55]) . With a gap this was done in pS] and without a gap 
similar arguments have been used in [14^ and in [15^ . We will not discuss this point 
further. Putting ([55]) and ([55)1 together, we reach the conclusion that Q ^ Q' 
except possibly on the kernel of Vq . Therefore 

for some < (5 < 1{o}(2^q)- If ker(2?Q) = {0} then of course 5 — Q. 

If ker(I?Q) 7^ {0} but has dimension > 2, it is a well-known fact that 5 = Q 
or (5 = 1{o}(2?q) that is, 5 must fill the last shell completely, see [21 13] and [23 
Prop. 3]. If dimker(I?Q) = 1, and Q + Pq is a projector, then, necessarily, (5 = 
or i5 = 1{o}(I'q) and we are done. If, however, Q + Pq is not a projector, then 
the argument does not work but, in this case, the energy does not change if we 
subtract 5 from Q, because the corresponding particle does not interact with itself, 
i.e., Q' = Q — 5 '\s a, minimizcr as well. Therefore we can redo the above argument 
with Q' instead of Q. But now, thanks to our definition p4p . we actually know 
that Q' + Pq is a projector. Hence, we deduce, as stated, that Q' = 1x{T>q>) — P^ 
with I = (-00,0) or I = (-00,0]. D 

Remark 6.1. Since Q must be a minimizcr for P^ . we deduce that 

lim fBDF(Qn) = £bDf{Q)- 

More precisely, we see that 

lim D{pq,^,pqJ^ D{pq,pq) 

n—>oo 

which implies that pq^^ — > pg strongly in C. Similarly we have 

hm tr p°|(g++ - Q--) = tr p°|(Q++ - Q-~) 

n— >-oo 

and 

lim // ^^^dxdy= lim // ^^dxdy. 

«^oo77jj2xR2 \x-y\ n^oo77jj2^]u2 \x - y\ 

By the reciprocal of Fatou's lemma for operators, this proves strong convergence of 
|X»°|^/^Q,t'^|27°|^/2 in ©1, and strong convergence of Qnix,y)\x ~ y\^^^^ in ©^ 
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Appendix A. Localization of massless kinetic energy 

In the following lemma, we provide an IMS-type formula for the massless rela- 
tivistic energy of a fermionic density matrix < 7 < 1, based on |33j . Results of 
the same form already exist in the literature, see in particular j391 Thm. 9]. There, 
a simple explicit formula for the operator \p\ — x\p\x ~ ^\'PV) is used. 

Lemma A.l. Let < 7 < 1 such that tr |p|7 < cx). Then for a partition of unity 
X^ + if' = 1, with smooth functions x,V one obtains in dimension d — 2 



2 

L4 



tr |p|7 > tr \p\xiX + tr \p\mV 

-citT\p\^f'{\\Vx\\h + \\yv\\hf'{\\yx\\ 
and in dimension d — 3 

tr \p\j > tr \p\x7X + tr b|?77?7 

- c (tr \p\^f' (II Vxlli. + II Vrylli.)'/' (|| Vx||ie 
where c is a universal constant. 
Proof. We use [33l Lemma A.l] which states that 



l^vUhf"- 



l^vUhf- 



\p\ - X\p\x - V\p\v > -; 

So, it remains to estimate the term 

1 



Vtdt^{\Vx\' + \Vvf) 



t + p'^ 



tr 



Vtdt- 



(|Vxl' + IVt^I' 



1 



Tl 



,0 t+p^ '' ' t+p^ 

We are going to decompose the integral into a small t and large t-part. For that 
reason we will use two different estimates for the integrand. Define 



Since < 7 < 1 we obtain 



(37) tr^-^/2^- 



r7<tr 



t+p^ t + p'^ 



f^{x)dx 



^^d/2 



(27r)''/2 
f^{x)dx 



dp 



{t+p^y 

dp 

(l+p2)2' 



on the other hand 
(38) tr 



71^/'7^7 = tr |p|-V2 1 j2 1 |^|-i/2|^|^ 
t + p'^ t + p^ t + p'^ t + p^ 



< 



1 



— r 



1 



tr |p|7 < 



1 



b|V 



1/2-^ 



tr |p|7. 



(i+p2)|p|l/2.' (i+p2)|p|l/2 

Recall that here / is seen as a multiplication operator on L'^{M.'^). In dimension 
d = 2 the Hardy-Littlewood-Sobolev inequality states 
2 



1 



\p\y 



1/2 •' 







= sup 




2(R2)- 


■>L2(I{2) 


ll'/'lli2(K2) 

= sup 


=1 






IIvL2{r2, 


=i 



{-A)-^/\f^) 



2 

L2(R2) 



f{x)ip{x)f{y)ip{y) 



< C sup ||/<y9||^4/3(R2) 

IIvIIl2(k2)=1 



\x-y\ 

-cii 



dxdy 



Il*(K2) 
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Combining these results we get in rf = 2 



< C (^V^||/||i.(K.) + C-i=||/||i.(R.) tr |p|7J 
<C||/|U2(K2)||/|U4(R.)(trb|7)i/^ 



after optimizing over e, which proves the Lemma in 2D. For d = 3 one proceeds 
similar with the corresponding Hardy-Littlewood-Sobolev inequality in 3D, leading 
to 

1 



w 



1/2 •' 



<c\\it 



L6(M3) 



This concludes the proof of Lemma lA.ll D 
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